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The source of human infection with Middle East respi¬ 
ratory syndrome coronavirus remains unknown. Molecular 
investigation indicated that bats in Saudi Arabia are infected 
with several alphacoronaviruses and betacoronaviruses. 
Virus from 1 bat showed 100% nucleotide identity to virus 
from the human index case-patient. Bats might play a role 
in human infection. 


S ince Middle East respiratory syndrome (MERS) was 
described in September 2012, over 90 cases have been 
reported worldwide, 70 from Saudi Arabia. The incidence 
of infection with the causative agent, a betacoronavirus 
(MERS CoV) (7), has not been determined; however, the 
mortality rate among those who received clinical care is 
*65% (2). Although instances of human-to-human trans¬ 
mission have been documented between case-patients and 
others in close contact (including hospital patients sharing 
rooms, family members, and medical personnel), the sourc¬ 
es of infection for most patients remain unknown. Because 
of sequence similarities between P-CoVs identified in bats 
and those of MERS CoV isolated from humans, a bat reser¬ 
voir has been posited (5-5). Although neither detection of 
MERS CoV in bats nor contact of human MERS patients 
with bats have been reported, a role for bats in human in¬ 
fection cannot be excluded because contact can be indirect 
(mediated through another animal vector or fomites). 
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The Study 

In October 2012 and April 2013, three agencies col¬ 
lected samples from bats in regions where MERS cases had 
been identified (Figure 1). The agencies are the Ministry of 
Health of Saudi Arabia, the Center for Infection and Im¬ 
munity of Columbia University, and EcoHealth Alliance. 

During the October investigation, the team inter¬ 
viewed the family of an index case-patient in Bisha and 
collected samples from bats <12 km from his home, in an 
abandoned date palm orchard, and <1 km from his place 
of employment, a hardware store that fronted a garden and 
date palm orchard. Although neither family members nor 
employees recalled seeing bats, the team observed roosting 
bats and guano in abandoned wells and ruins within 12 km 
of his home and insectivorous bats at dusk in the garden 
behind his store. Over 3 weeks, 96 bats representing 7 spe¬ 
cies {Rhinopoma hardwickii, Rhinopoma microphyllum, 
Taphozous perforatus, Pipistrellus kuhlii, Eptesicus bot- 
tae, Eidolon helvum, and Rosettus aegyptiacus) were cap¬ 
tured in mist nets and harp traps, then released after visual 
speciation and collection of morphometric measurements; 
wing punch biopsy samples; blood; throat swab samples; 
and rectal swab samples or fecal pellets. Samples were col¬ 
lected into viral transport medium or lysis buffer. 

During the 3-week April investigation, fecal samples 
were collected on tarps laid out at bat roosting sites in and 
around Bisha, Unaizah, and Riyadh. Representative animals 
at each roosting site were captured, identified morphologi¬ 
cally, and released after wing punch biopsy samples were 
collected for speciation by DNA analysis. Samples were 
collected into cryovials. 
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Figure 1. Bat sampling sites and locations of home and workplace 
of index case-patient with Middle East respiratory syndrome, Bisha, 
Saudi Arabia. 
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DISPATCHES 




Table 1. PCRs and primers used in CoV detection* 



PCRs (reference) 

Primers, 5’->3’ 

Nested fragment size, 
region (primer locations on 
the reference genome)f 

Type of 
CoV (no.) 

Nested pan-CoV-l 
(6) 

PLQ-F1, CGTTGGIACWAAYBTVCCWYTICARBTRGG 
PLQ-R1, GGTCATKATAGCRTCAVMASWWGCNACATG 
PLQ-F2, GGCWCCWCCHGGNGARCAATT 

PLQ-R2, GGWAWCCCCAYTGYTGWAYRTC 

*400 nt, RdRp 
(18310-187450) 

a-CoV (8), 
P-CoV (1) 

Nested pan-CoV-ll 
(7) 

WT-C0V-F1, GGTTGGGAYTAYCCHAARTGTGA 
WT-COV-R1, C CATC AT CAS W YRAAT CAT CATA 
WT-COV-F2, GAYTAYCCHAARTGTGAYAGAGC 
WT-COV-F3, GAYTAYCCHAARTGTGAUMGWGC 

*430 nt, RdRp 
(15260-15700) 

a-CoV (5), 
P-CoV (2) 

Hemi-nested RdRp- 
sequence assay (9) 

EMC-SeqRdRP-Rev, GCATWGCNCWGTCACACTTAGG 
EMC-SeqRdRP-Fwd, TGCTATWAGTGCTAAGAATAGRGC 
EMC-SeqRdRP-Rnest, CACTTAGGRTARTCCCAWCCCA 

*230 nt, RdRp 
(15048-15290) 

a-CoV (2), 
P-CoV (1) 

Hemi-nested N- 
sequence assay (9) 

EMC-SeqN-Fwd, CCTTCGGTACAGTGGAGCCA 
EMC-SeqN-Rev, GATGGGGTTGCCAAACACAAAC 
EMC-SeqN-Fnest, TGACCCAAAGAATCCCAACTAC 

*280 nt,N seq 
(29,549-29,860) 

NA 

Nested Cll-pan- 
CoV-lll 

NM-CoV-2F1, ACWGTTCARGGICCWCCIGG 
NM-CoV-2F2, GTTCARGGGCCWCCGGGNAC 
NM-CoV-2R1, GGCAGCTGWGCWGGRTCICCNACRTA 
NM-CoV-2R2, AGCTGWGCWGGRTCGCCIACRTANAC 

*355 nt, helicase 
(17,060-17,410) 

P-CoV (2) 

Nested CII-MERS- 
RdRp 

NM-HCOV-F1, GTGCTAAGAATAGAGCTCGCACT 
NM-HCOV-F2, AGAGCTCGCACTGTTGCAGGC 
NM-HCOV-F2, AGAGCTCGCACTGTTGCAGGC 
NM-HCOV-R1, ACCCAT AAGAT GCGG ATT AT CAAC 
NM-HCOV-R2, TGCGGATTATCAACATCTTTGTAC 

*190 nt, RdRp 
(15068-15249) 

P-CoV (1, 
MERS 
CoV) 

Hemi-nested Cll- 
MERS N sequence 

NM-NSeq-F-1, ACTTCCTTCGGTACAGTGGAGC 
NM-NSeq-R-1, GGCACTGTTCACTTGCAATC 
NM-NSeq-R-2, GGAGGTTCAGACATTTGGTCT 

*170 nt, N seq 
(29545-29713) 

NA 

upE and ORFIb 
real-time assays (8) 

upE-Fwd: GCAACGCGCGATTCAGTT 
upE-Prb: FAM-CTCTTCACATAATCGCCCCGAGCTCG-TAMRA 
upE-Rev: GCCTCTACACGGGACCCATA 

ORFIb-Fwd: TTCGATGTTGAGGGTGCTCAT 

ORF1 b-Prb: FAM-CCCGTAATGCATGTGGCACCAATGT-TAMRA 
ORF1 b-Rev: TCACACCAGTTGAAAATCCTAATTG 

Upstream of E gene and 
ORF 1b 

NA 

*CoV, coronavirus; MERS, Middle East respiratory syndrome; RdRp, RNA-dependent RNA polymerase; NA, not applicable; ORF, open reading frame. 
tPrimer locations are based on human (3-CoV 2c EMC/2012, complete genome (GenBank accession no. JX869059). 


All samples were stored in liquid nitrogen and con¬ 
veyed to Riyadh for storage at -80°C before being trans¬ 
ported to Columbia University in New York in dry nitro¬ 
gen. The October 2012 shipment was inadvertently opened 
at customs in the United States and sat at room temperature 
for 48 hours before transfer to Columbia University; at ar¬ 
rival, all samples had thawed. The April 2013 samples ar¬ 
rived intact. 

Total nucleic acid was extracted from samples by us¬ 
ing the NucliSENS easyMAG system (bioMerieux, Dur¬ 
ham, NC, USA) and subjected to 8 PCRs with primers and 
protocols designed to amplify regions within the helicase, 
RNA-dependent RNA polymerase (RdRp), and nucleo- 
capsid or envelope proteins of CoVs (6-9). Products were 
sequenced and analyzed for similarity to GenBank data¬ 
base entries by using the BLASTn and BLASTx programs 
(www.ncbi.nlm.nih.gov/blast/Blast.cgi). Primer sequences 
are shown in Table 1. The identity of bat species yielding 
specific viral products was determined by amplifying and 
sequencing a fragment of the cytochrome B gene (10). All 
visual classifications of species were confirmed except for 


that of T. perforatus bats. There is no reference sequence 
for T. perforatus bats in GenBank. However, because the 
closest reference sequence was from T. nudiventris bats, at 
84% identity we presume that the product represents bona 
fide T. perforatus bat cytochrome B gene sequence. Repre¬ 
sentative cytochrome B sequences have been uploaded to 
GenBank (accession nos. KF498635-KF498641). 

Table 1 indicates the CoV genera identified by using 
individual primer sets. As anticipated, pan-CoV assays de¬ 
tected a- and P-CoVs. One assay specific for MERS CoV 
(9) also detected a-CoVs. This finding reinforces the need 
for sequence confirmation of PCR products. Table 2 indi¬ 
cates the CoV species identified with respect to location, 
sample type, and bat species. CoV sequences were ampli¬ 
fied from rectal swab samples or fecal pellets and from 
roost feces but not from serum, throat swab samples, or 
urine. a-CoV sequences were amplified more frequently 
than P-CoV sequences (223 vs. 4). Whereas a- and P-CoV 
sequences were amplified from CoVs from T. perforatus, 
E. helvum, and R. hardwickii bats, only alpha sequences 
were amplified from CoVs from P. kuhlii bat samples. 


1820 


Emerging Infectious Diseases 'www.cdc.gov/eid • Vol. 19, No. 11, November 2013 









MERS CoV in Bats, Saudi Arabia 

Table 2. CoVs detected in bats, Saudi Arabia* 








No. samples tested (no. positive) 

Total no. 


Bat family, genus, 
species 

No. 

Location bats 

Throat 

swab 

Fecal 

pellets 

Roost 

Urine Serum feces 

samples, 
n = 1,003 

Total no. positive samples, n = 
227 (closely related CoVs)f 


October 2012 


Emballonuridae 


Taphozous perforatus 

Bisha ruins 

29 

29 (0) 

25 (2) 

8(0) 

22 

10 (1) 

94 

1 p-Cov (1 MERS novel CoV) 
and 2 a-CoVs (1 bovine 
respiratory CoV, 1 Kenya bat 
CoV BtKY86) 

Pteropodidae 

Eidolon helvum 

Bisha town 
center 

25 

25 (0) 

25 (5) 

13(0) 

19 

NA 

82 

1 p-Cov (1 Eidolon bat CoV- 
HKU1) and 4 a-CoVs (4 Kenya 
bat CoV BtKY86) 

Rousettus 

aegyptiacus 

Bisha town 
center 

3 

3(0) 

3(0) 

1 (0) 

2 

NA 

9 

NA 

Rhinopomatidae 

Rhinopoma 

hardwickii 

Naqi and 
Old Naqi 

36 

36 (0) 

35 (0) 

4(0) 

NA 

15(0) 

90 

NA 

Rh. microphyllum 

Old Naqi 

1 

1 (0) 

1 (0) 

NA 

NA 

NA 

2 

NA 

Vespertilionidae 

Eptesicus bottae 

Bisha ruins 

1 

1 (0) 

1 (0) 

1 (0) 

NA 

32 (0) 

35 

NA 

Pipistrellus kuhlii 

Bisha ruins 

1 

1(0) 

1 (0) 

NA 

NA 

NA 

2 

NA 


April 2013 

Rhinopomatidae 


Rh. hardwickii 

Greater 
Bisha area 

NA 

NA 

NA 

NA 

NA 

209 

(93) 

209 

2 p-Covs (2 canine respiratory 
CoVs) and 91 a-CoVs (5 canine 
CoVs, 2 Miniopterus bat CoVs, 
84 Chaerephon bat CoV) 

T. perforatus 
Vespertilionidae 

Bisha ruins 

NA 

NA 

NA 

NA 

NA 

203 (0) 

203 

NA 

P. kuhlii 

Greater 

Unaizah 

area 

9 

9(0) 

NA 

NA 

NA 

263 

(126) 

277 

126 a-CoVs (69 alphaCoV 
P.kuh-Spain, 3 canine CoVs, 

37 bat CoV P.pyg/Germany, 1 
human CoV NL63, 2 Rousettus 
batCoVHKUIO, 11 porcine 
epidemic diarrhea virus, 2 
Cardioderma bat CoVs, 1 
Hipposideros bat CoV HKU10) 


Greater 

Riyadh 

area 

5 

5(0) 

NA 

NA 

NA 

NA 

NA 

NA 


*CoV, coronavirus; MERS, Middle East respiratory syndrome; NA, not applicable. 
fBased on BLASTn (www.ncbi.nlm.nih.gov/blast/Blast.cgi). 


CoV sequences were amplified from 220 of 732 roost 
feces samples and 7 of 91 rectal swab samples or fecal pel¬ 
lets. A product obtained by PCR amplification of nucleic 
acid from a fecal pellet of a T. perforatus bat captured in 
October 2012 in Bisha showed 100% nt identity to the hu¬ 
man P-CoV 2c EMC/2012 cloned from the index case-pa¬ 
tient in Bisha. A phylogenetic analysis of CoVs obtained 
in this study is shown in Figure 2. CoV sequences have 
been uploaded in GenBank (accession nos. KF493884- 
KF493888). 

Conclusions 

A wide range of CoV species are circulating among 
bats in Saudi Arabia. Although the prevalence of CoVs 
was high (=28% of fecal samples), MERS CoV was found 
in only 1 bat. A 3.5% MERS CoV infection rate (n = 29; 
95% Cl 0-20%) in T. perforatus bats is low compared 


with that for severe acute respiratory syndrome-like CoV 
in rhinolophid bats in China (10%—12.5%) but consistent 
with CoV prevalence among bats in Mexico (4). Further¬ 
more, the sensitivity for viral nucleic acid detection in 
samples collected in October 2012 was probably reduced 
because of failure in cold chain transport. Whereas 219 
(32%) of 675 of fecal pellets collected in April revealed a 
CoV sequence by PCR, only 8 (5%) of 148 of rectal swab 
samples or fecal pellets collected in October were positive 
by the same assays. We were unable to recover additional 
sequences beyond the 190-nt RdRp fragment represented 
in Figure 2 but are confident in the fidelity of the finding. 
First, although RdRp is a conserved portion of the CoV 
genome, there is no precedent for 100% identity of a bat 
sequence with a human MERS CoV sequence. Second, 
when this work began we did not have cultured MERS 
CoV, human MERS samples, or MERS CoV cDNA in the 
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DISPATCHES 


KSA-287- Taphozous perforates bat 
JX869059-Human betacoronavims 2c EMC/2012 
KC164505-Betacoronavirus England 1 
KC667074-Human betacoronavims 2c England-Qatar/2012 
KF186567-Middle East respiratory syndrome coronavirus-AI-Hasa 1 2013 
KF192507-Middle East respiratory syndrome coronavims 
KC776174-Human betacoronavims 2c Jordan-N3/2012 
KC243390-Bat coronavirus BtCoV/8-724/Pip pyg/ROU/2009 
KC243392-Bat coronavirus BtCoV/UKR-G 17/Pip nat/UKR/2011 
KC243391-Bat coronavims BtCoV/8-691/Pip nat/ROU/2009 
KC522088- Pipistrellus bat coronaurus HKU5 isolate 33R 
KC522088- Pipistrellus bat coronaurus HKU5 isolate 32R 

_j- KC522044-7y/oncyfens bat coronavims HKU4 isolate 10R 

L EF065505-Bat coronawrus HKU4-1 
-HQ728482-E/do/on bat coronawrus/Kenya/KY24/2006 


MERS CoV 


99 


98 


_99j AB514460-Human coronavims HKU1-ORF1b-PHL07-1656 
■ DQ415911-Human coronavims HKU1-N15 genotype B 
JN129834-Human coronavirus OC43-HK04-01 
AB514461-Human coronavims OC43-ORF1b-PHL07-1581 


Ir FJ938065-Boune respiratory coronavims-AH187 
> KSA-004- Rhinopoma hardwickii bat 
■ JX860640-Canine respiratory coronavims-K37 
|— FJ588686-Bat SARS Co-Rs672/2006 
4- GQ153547-Bat SARS coronavirus-HKU3-12 
99 L GQ153540 Bat SARS coronavirus HKU3-5 


JX503060-Human coronavirus 229E isolate 0349 

86j KC175339.11 Canine coronawms strain 171 complete genome 
AY994055.1| Feline infectious peritonitis virus 
DQ811787-PRCV-ISU-1 

HQ012368-Feline coronavirus-UU18 
FJ938061-Feline coronavims-UU3 


86. KC1753: 
92 f AY9940J 
83 I— DQ81 

73 I_ I - 

jT 1 FJ93 


- HM245925-Mink coronavirus Strain-W01127 

- DQ249235-Bat coror.a^rus-HKU2-1 
HQ728485-M/mopferus bat coronavirus/Kenya/K Y33/2006 

- EU420138-Bat coronavims 1A-AFCD62 
100 | KSA-002-R/)/nopoma hardwickii bat 

• K SA-001 -Plpls trellus-kuhlii bat 
HQ728486-C/7aerephon bat coronavims/Kenya/KY22/2006 

- KS A-OOZ-Rhinopoma hardwickii bat 

I-AF353511 PEDV strain CV777 

I-HQ728480-Carc//oderma bat coronavirus/Kenya/KY43/2006 

JQ989266-H/ppos/denos bat coronavims HKU10 isolate TT3A 

-EU769557-Bat coronavims Trinidad/IFY2BA/2007 

I-EU420139-8at coronavims HKU8 strain-AFCD77 

•— AY567487-Human coronavims NL63 


Figure 2. Phylogenetic tree 
showing genetic relatedness 
between coronaviruses 

identified in bat samples from 
Saudi Arabia (boldface), 
MERS coronaviruses, and 
other published coronavirus 
sequences available in 
GenBank. The maximum- 
likelihood tree of partial RNA- 
dependent RNA polymerase 
gene (nt position 15068-15249 
of GenBank accession no. 
JX869059) was constructed 
by using the Tamura-Nei 
model with discrete gamma 
rate differences among sites 
as implemented in MEGA 5.2 
(www.megasoftware.net). Each 
branch shows the GenBank 
accession number followed by a 
brief description of the sequence 
used. Scale bar indicates 
nucleotide substitutions per site. 
MERS, Middle East respiratory 
syndrome; CoV, coronavirus; 
SARS, severe acute respiratory 
syndrome; KSA, Kingdom of 
Saudi Arabia. 


laboratory at Columbia University where samples were 
removed directly from the tubes in which they were col¬ 
lected in the field for nucleic acid extraction, PCR, and 
sequence analysis. Third, the only MERS-positive signal 
was obtained in PCR analysis of the T. perforatus bat cap¬ 
tured in Bisha near the home and workplace of the MERS 
index case-patient used to generate the human P-CoV 2c 
EMC/2012 sequence. 

Bats are reservoirs of several viruses that can cause hu¬ 
man disease, including rabies, Hendra, Nipah, Marburg, se¬ 
vere acute respiratory syndrome CoV, and probably Ebola 
viruses (11-14). Cross-species transmission from bats to 
humans can be direct, through contact with infected bats or 
their excreta, or facilitated by intermediate hosts (15). Bat 
CoVs are typically host specific; however, MERS-related 
Co Vs have reportedly been found in many bat families, in¬ 
cluding Vespertillionidae, Molosidae, Nyteridae, and now 
Emballonuridae (sheath-tailed bats) in Africa, the Ameri¬ 
cas, Asia, and Europe. We sampled only a small sample 
of bats in Saudi Arabia. Nonetheless, given the rarity of 
MERS CoV sequences detected by our survey and the 
broad distribution of MERS cases throughout the Middle 


East, we speculate that there are probably other hosts. Fu¬ 
ture work should investigate additional bat and other wild¬ 
life species and domestic animals for CoV infection and 
potential linkage to human disease. 
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